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Chapter  1 
INTRODUCTION 


The  turbulent  nature  of  the  atmosphere  produces  random  fluc¬ 
tuations  in  Its  refractive  Index  fields.  The  refractive  Index 
perturbations  cause  .  degradation  of  the  Initial  phase  coherence 
of  optical  beam,  vhlch  propagate  through  the  turbulent  atmospheric 
medium.  This  phase  coherence  degradation  can  limit  the  performance 
of  precision  pointing  and  tracking  systems  and  high  resolution 
optical  systems  vhlch  must  operate  over  an  atmospheric  propagation 
path.  Statistical  data  describing  the  phase  coherence  variations 
induced  by  the  atmosphere  Is  therefore  repaired  for  performance 

limit  evaluation  and  error  analysis  for  the  previously  described 
systems. 

A  convenient  and  useful  measurement  of  the  phase  degradation 

is  the  measurement  of  the  temporal  fluctuations  of  the  optical  phase 

difference  between  two  points  located  at  the  entrance  plane  of  a 

receiver.  The  experimental  measurements  were  performed  by  an 

optical  interferometer  which  measured  the  phase  difference  between 

the  two  sampling  apertures.  The  experiment  consisted  of  measuring 

Phase  difference  fluctuations  for  various  propagation  conditions 

The  phase  difference  data  was  then  computer  processed  by  fast 

tier  transform  techniques  to  yield  experimental  power  spectral 
density  estimates. 

The  characterization  of  the  observed  phase  difference  fluc- 
'  -  -1'ns  in  terms  of  their  power  spectral  density  yields  Information 
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directly  applicable  to  design  and  error  analysis  studies  for  optical 
wavelength  systems  operating  in  the  atmospheric  medium.  Graphs  of 
the  experimental  phase  difference  power  spectral  density  are  compared 
to  graphs  of  theoretical  power  spectral  density  functions  derived 
from  statistical  models  of  the  atmospheric  refractive  index  fields. 

A  recently  proposed  model  which  agrees  with  the  experimental  spectral 
data  at  all  frequencies  is  described. 

As  an  important  backup  to  the  phase  measurements,  microtempera¬ 
ture  spectra  were  measured  during  the  course  of  optical  missions. 
Temperature  fluctuations  govern  refractive-index  fluctuations  in  the 
visible  to  middle  infrared  wavelengths.  Thus,  any  discrepancies  between 
theoretical  and  experimental  phase  spectra  should  also  be  seen  in  the 
microtemperature  spectra.  In  fact,  all  measured  spectra  demonstrate 
more  power  in  the  low  frequencies  than  is  indicated  by  theory. 


Chapter  2 

THEORETICAL  BACKGROUND 

2.1  Introduction 

Optical  beams  propagating  through  the  atmosphere  undergo 
degradation  of  their  phase  and  amplitude  characteristics  The 
degradation  arises  from  the  fact  that  the  Index  of  refraction 
fields  of  the  atmosphere  are  statistical  random  functions  of  their 
spatial  and  temporal  coordinates.  In  recent  years  major  emphasis 
has  been  placed  on  understanding  the  effect  of  the  random  refractive 
Index  fields  on  beam  propagation  for  laser  wavelengths.  Particular 
attention  has  been  paid  to  phase  characteristics  because  of  the 
coherent  nature  of  laser  light  and  the  need  to  maintain  coherence 
for  a  number  of  laser  applications. 

Previous  experiments  have  been  conducted  to  determine  spatial 
and  temporal  characteristics  of  laser  beams  operating  at  visible 
wavelengths.  The  experiments  reported  here  were  aimed  at  deter¬ 
mining  the  optical  phase  difference  at  two  points  across  a 
spherical  wave,  10.6-mlcron  wavelength  laser  beam  whose  Initial 
phase  coherence  was  destroyed  by  propagation  through  the  atmosphere. 
Theoretical  expressions  for  the  phase  difference  power  spectrum 
have  been  derived  by  several  Investigators.  The  differences  In 
these  various  expressions  for  the  spectrum  are  due  to  the  different 
assumptions  used  In  their  derivation.  This  section  will  present 
several  theoretical  expressions  for  the  phase  difference  power 
spectrum  and  their  most  important  characteristics,  underlying 
assumptions,  and  differences. 
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Several  basic  assumptions  about  the  characteristics  of  the 
refractive  index  fields  of  the  atmosphere  are  common  to  all  the 
derivations  and  resulting  expressions.  These  assumptions  will 
be  explained  before  presenting  the  different  expressions.  Two 
assumptions  which  simplify  the  theory  are  that  the  atmospheric 
refractive  index  fields  are  three-dimensionally  isotropic  and 
homogeneous .  The  assumption  of  isotropy  simplifies  the  three- 
dimensional  spatial  correlation  functions  of  the  refractive 
index  field  to  functions  of  single  variable  scalar  arguments. 

The  assumption  of  statistical  homogeneity  of  the  atmospheric 
index  fields  permits  writing  the  two-point  spatial  correlation 
functions  as  functions  of  a  single  variable,  the  separation  of 
the  two  observation  points.  Also  used  for  the  derivation  cf  all 
the  following  phase  difference  power  spectrum  equations  i3  the 
assumption  that  temporal  changes  in  the  atmospheric  refractive 
index  fields  at  any  point  can  be  explained  by  use  of  Taylor's 
hypothesis. 

Taylor's  hypothesis  proposes  that  the  large  scale  mean 
atmospheric  motions,  such  as  those  produced  by  wind,  are  the 
dominate  transport  factor  for  moving  "frozen  in"  spatial  varia¬ 
tions  of  the  refractive  index  field  past  the  observation  point. 
Through  the  use  of  this  model,  knowledge  of  the  spatial  statistics 
of  the  index  field  can  be  transformed  to  temporal  statistics  by 
use  of  the  wind  speed. 
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2.2  Theoretical  Phase  Difference  Power  Spectra 


Tatarski^  published  the  first  derivation  of  an  equatic 


for  the  phase  difference  power  spectrum  of  a  light  beam  propa¬ 
gating  through  a  turbulent  atmosphere.  This  expression  was  for 
the  special  case  of  plane  wave  propagation  and  therefore  does 
not  directly  apply  to  this  experiment  which  was  conducted  with 


spherical  waves. 

In  1971  Clifford^  published  the  theory  for  the  spherical 
wave  phase  difference  power  spectrum.  Clifford's  resulting 
expression  for  W(f ) ,  the  phase  difference  power  spectrum,  is 

W(f)  "  [2*0]  (*033)(2irp)^k2LCn  “2r^"]  **  (2_1) 


with  the  restrictions  that  f>(vL/2lTL0)  and  S-0«XL«L0 
where  f  -  frequency  in  hertz 

k  -  2  r/X,  the  wavenumber  of  the  light  beam 
X  *■  wavelength 
L  ■  propagation  path  length 

C2  *  the  optical  refractive  index  structure  constant 
n 

p  -  separation  of  the  two  phase  difference  observation  apertures 
v.  ■  wind  velocity  perpendicular  to  the  propagation  path 

•  k 

L0  ■  size  of  the  outer  scale  of  turbulence 

JL  ■  size  of  the  inner  scale  of  turbulence 
0 

X  -  (2TTfp/Vj_) 

1.0  for  p«v/u 
2.0  for  p»y^XL’ 


Graphs  of  the  above  equation  for  several  typical  wind  velocities 
are  shown  in  Fig.  2.1.  The  curves  have  a  -8/3  power  slope  at  higher 
frequencies,  f>  VTv  /2-np,  and  a  -2/3  slope  at  lower  frequencies 
where  v^/2irLo  <f< v^/2np.  The  curves  do  not  extend  below  the 
frequency  f  ■  v^/2ttL0  because  that  is  the  lower  limit  of  validity 
of  the  equation. 

Clifford's  derivation  of  the  phase  difference  power  spectrum 
Included  the  three  previously  explained  suppositions  of  three- 
dimensional  Isotropy,  statistical  homogeneity,  end  applicability 
of  Taylor's  hypothesis  for  the  atmospheric  refractive  index  fields. 
Clifford  also  assumed  that  the  vector,  p,  which  connects  the  two 
phase  difference  observation  apertures  and  the  vector, v^,  which 
describes  the  horizontal  wind  velocity  are  coplanar  in  a  horizontal 
plane  and  are  also  parallel.  Implicit  in  Clifford's  derivation  is 
the  applicability  of  the  Kolmogorov  refractive  index  spectrum 

•  (K)  -  .033  C2  K_ll/3  (2-2) 

n  n 

in  which  the  variable  K  ■  2 is  the  spatial  frequency  proportional 
to  scale  size  £..  This  spectrum  is  valid  only  for  refractive  index 
perturbations  of  scale  sizes  in  the  range  £0<£<L0  where  and  Lo 
are  the  inner  and  outer  scales  of  turbulence  in  the  atmosphere. 

The  use  of  this  specific  refractlvity  spectrum  is  the  source  of 

the  two  restrictions,  f>v^/2TTL0  and  io<<V^L  «Lo,  on  Clifford's  final 

equation. 

Collins^)  has  presented  two  different  expressions  for  the 
phase  difference  power  spectrum.  The  first  expression  extends 


the  region  of  validity  of  Clifford's  spectrum  to  frequencies  less 
than  v1/2ttLo  and  is  for  the  case  where  the  vectors  p  and  v^are 
horizontally  coplanar  and  parallel.  The  second  expression  is  also 
valid  for  frequencies  less  than  v,/2ttLa  but  is  for  the  case  for 
which  vectors  p  and  v^  are  not  horizontally  coplanar  and  parallel. 

Collins'  phase  difference  power  spectrum  equation  for  the 
case  in  which  vectors  p  and  v^  are  horizontally  coplanar  and  parallel 

is 

«(f)  -  [A  (.0330  (2  Trp)  *  k«u£  »r  l**  +  -f*-]  (2'3) 

with  the  restriction  that  2-o<<  where  R  ■  Ap/L0  ,  A  »  1.071,  and 

all  other  quantities  are  as  defined  for  (2-1).  Fig.  2.2  shows 
this  equation  plotted  for  a  variety  of  typical  values  of  ^ ,  the 
mean  wind  speed  perpendicular  to  the  propagation  path  and  horizontally 
coplanar  with  and  parallel  to  vector  p.  The  spectra  in  Fig.  2.  2 
are  for  a  small  observation  aperture  separation,  p,  and  have  a  -8/3 
power  dependency  in  the  high  frequency  region  where  v^/2*P , 

a  -2/3  power  slope  in  the  mid-frequency  region,  and  a  +2  slope  for 
the  low  frequency  region  where  f<  v^/2ttL0. 

Collins  has  also  derived  an  equation  for  the  phase  difference 
power  spectrum  for  the  case  in  which  the  mean  horizontal  wind  velocity 
vector,  v^,  and  the  observation  aperture  separation  vector,  p,  are 
not  parallel  and  coplanar  in  a  horizontal  plfne.  The  geometry  for 
the  vectors  p  and  v^  ,  as  shown  in  Fig.  2.3,  is  such  that  when  viewed 
in  a  plane  normal  to  the  propagation  path  they  are  separated  by  an 
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APERTURE  2 


orientate ^  angle,  0.  The  resulting  power  spectral  density  equation 
has  two  general  forms;  the  first  is  valid  for  low  and  mid-range 
frequencies  and  tb&  second  for  nigh  frequencies.  The  two  forms  and 
their  regions  of  validity  are 

w(f)  “|*o]  (*033)  (2irp)  ^l^LC^v^1  jx2+R^j”^3  j*  c°a.  (-9)  +  (X..  +*  iat1L.L9).j(2-4) 

valid  for  |xcos(0)  +  (X2+R2)  ^2sin(0)  J<V^* 
and 


W(f) 


(.(mHZTrp)8^2!^2^1  [x2+r2]"^ 


(2-5) 


valid  for  |Xcos(0)  +  (X2+R2)  ^sin(0)  J  > 


The  above  expressions  properly  combined  to  extend  over  both  regions 
of  validity  are  shown  in  Fig.  2.4  for  small  separation,  p,  the 
case  of  interest  for  this  experiment.  The  high  and  medium  frequency 
regions  are  independent  of  the  orientation  angle  and  exhibit  the 
previously  seen  -8/3  and  -2/3  power  characteristics,  respectively. 
However,  in  the  low  frequency  region,  f<v^/2irL0,  the  curves  flatten 
out  at  a  value  proportional  to  the  quantity 


(.033)  k2LC2 


(2irp)  ^v^1  (0,25  Fl  sin20) . 


The  basic  assumptions  used  by  Collins  to  derive  the  theoretical 
phase  difference  spectra  for  both  cases  are  identical.  The  first 
three  suppositions,  also  used  by  Clifford,  are  three-dimensional 
isotropy,  homogeneity,  and  the  applicability  of  Taylor's  hypothesis 
for  the  atmospheric  refractive  index  fields.  For  both  of  Collins' 
derivations  the  Kolmogorov  refractive  index  spectrum  used  by 
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Clifford,  was  changed  to  the  modified  von  Karman  spectrum, 

0  (K)  =  (.033)0’  f(A/L0)  2+K2]  “ 1  1/6 .  (2-6) 

n  n  l  J 

The  parameter  K  »  2n/£  is  the  refractive  index  spatial  frequency 
corresponding  to  scale  size  A.  In  the  spatial  domain  this  refrac- 
tivity  spectrum  is  theoretically  valid  for  scale  sizes  both  within 
the  inertial  subrange  of  turbulence,  £,o<£<L0,  an^  f°r  scale 
sizes  larger  than  the  outer  scale  of  turbulence,  £>Lo .  This  exten¬ 
sion  of  the  refractive  index  spatial  spectrum  to  include  scale  sizes 
greater  than  the  outer  scale  was  responsible  for  extending  the  temporal 
phase  difference  spectrum  to  the  low-frequency  region  defined  as 
f <v^/2nL0 . 

Greenwood has  performed  an  analysis  of  the  phase  difference 
power  spectrum  which  is  an  extension  of  Collins'  theory  for  a  nonzero 
orientation  angle.  Greenwood's  analysis  assumes  that  the  orientation 
angle,  0,  is  a  random  variable  which  is  dependent  upon  the  horizontal 
wind  direction,  the  vertical  wind  velocity  variance,  and  the  small 
offset  angle  of  the  observation  apertures  from  the  horizontal  position. 
The  final  result  is  a  mean  phase  difference  power  spectrum,  <W(f)>, 
averaged  over  the  random  variable  0,  the  orientation  angle.  The 
equations  for  the  spectrum  are 

<W(f)>  -jl.oj  (.033)(2TTP)^k2LC’v[1[x2+R]'‘y3  |^j  l+ffrj 
\alid  for  |Xcos(6)  +  (X2+R2)  ^sin(6)  | « V&” 
and 

<W(f )>  -  ^'0  (•°33)(2^P)^k2LC2v[I[x2+R2]_v/j  (2-8) 

valid  for  |Xcos(d)  +(X2+R2)  ^sin(9)  |  »' 
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where  r  -  <02>  -  <0>2  +  <  6-<e>  2>* 

The  brackets  indicate  ensemble  averages  over  the  random  variable  0 

2 

Replacing  the  X  / 6  term  in  the  first  expression  above  by 

tl-(sinX) /X]  results  in  the  expression 
<W(f)>  -  (.033)  (2irp)  *^k2LC*v^ 1  [X2+R2]-1^' 

which  is  a  generalization  of  (2-3)  for  all  frequencies. 

The  parameters  in  the  above  liquation  which  have  the  greatest 
effect  on  the  shape  of  the  spectrum  are:  r,  the  second  moment  of 
the  orientation  angle;  L0,  the  turbulence  outer  scale  size  which 
determines  R;  and  v^ ,  the  mean  horizontal  wind  speed  perpendicular 
to  the  propagation  path.  Fig.  2.5  is  a  graph  of  (2-9)  for  values 
of  1,  2.5,  5,  7.5,  and  10  meters  per  second  for  parameter  v^ .  For 
larger  values  of  the  shift  to  the  right  of  the  upper  and  lower 
frequency  breakpoints  of  the  spectrum,  f  -rfb  and  f  -  A^/2ttL0, 

respectively,  is  very  clearly  shown.  Fig.  2.6  is  a  graph  of  the 
spectrum  for  various  typical  values  of  parameter  r,  for  an  outer 
scale  of  one  meter,  and  for  a  normal  mean  wind  speed  of  one  meter 
per  second.  Fig.  2.7  is  another  graph  for  the  same  parameters  as 
Fig.  2.6  with  the  exception  of  outer  scale  size  which  is  two  meters. 
These  two  preceding  graphs  clearly  show  the  effect  of  parameter  r 
on  the  shape  of  the  spectrum  at  low  frequencies.  Fig.  2.8  shows 
that  the  effect  of  an  increasing  outer  scale  size  on  the  low 
frequency  portion  of  the  spectrum  is  to  lower  the  low  frequency 
spectral  breakpoint  of  f  -  Avl/2ttL0  and  raise  the  low  frequency 
level  of  the  curve. 
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PRBQUENCY  (HERTZ) 

Theoretical  phase  difference  spectra  derived  by  Greenwood  for 


(ziHaH/ySNViaro) 
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In  general  the  equation  for  the  phase  difference  power  spectrum 
developed  by  Greenwood  has  the  -8/3  and  approximate  -2/3  slopes  at 
high  and  medium  frequencies  as  exhibited  by  the  previous  spectra 
derived  by  Clifford  and  Collins.  However,  unlike  the  previous  curves 
for  small  orientation  angles,  these  curves  show  a  low-frequency  flat¬ 
tening  trend  at  a  level  within  one  decade  of  the  medium  frequency 
spectral  peak.  The  actual  level  of  the  theoretical  spectra  at  low 
frequencies  is 


(.033)k2Ix£(2irp)**v£l(l/4)Ra. 

As  previously  stated  Greenwood's  expression  is  a  modification 
of  Collins'  work,  equations  2-4  and  2-5,  for  a  nonzero  orientation 
angle,  0.  The  assumptions  employed  by  Greenwood  include  all  those 
used  by  Collins  plus  several  additional  assumptions  explained  below. 
The  orientation  angle,  8,  is  assumed  to  be  a  Gaussian  distributed 
random  variable,  thereby  simplifying  the  determination  of  mean  values 
for  cos20  and  sin28  according  to  the  equations 


<co82(9)>  ■  ” 
<sin2(0)>  -  i 


1  +  e~za  cos(2n) 


1  -  e  20  cos(2h) 


(2-10) 

(2-11) 


where  H  *  the  mean  value  of  8 
o  ■  the  variance  of  8. 

For  this  experiment  we  may  assume  that  the  mean  value  of  0  is  very 
small  and  that  the  variance  of  0  is  also  fairly  small  so  that  the 
above  equations  may  be  approximated  by 


<cos2(0)>=  1 
<sin2 (0) >  =  <02>. 
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(2-12) 

(2-13) 


Fig.  2.9  shows  the  geometry  describing  the  resulting  situation. 
Restricting  the  solution  to  small  values  of  9  such  that  Q  ~  tan6 


and  assuming  a  small  variance  for  the  component  of  horizontal 
wind  speed  perpendicular  to  the  propagation  path  the  resulting 
approximation  is 


or 


<e2>  = 


<Vvert> 

<vx>2 


(2- 


<02>  s<Vvert>  _ 

<V>2cOS'  (<+>) 

where  Vvert  »  the  vertical  wind  speed  component 

v^«=  the  horizontal  wind  speed  component  perpendicular 
to  the  propagation  path 
V  ■  the  horizontal  wind  speed 
0  *  the  angle  separating  <V>  and  <v.  >. 


(2- 


From  boundary  layer  turbulence  theory  the  expression  for  estimating 

(5,6) 

the  second  moment  of  the  vertical  wind  velocity  is  given  as 


<Vv«rt>=  1*75 


[■4  <V>  1 

ln(Z/Zoy 


(2- 


or 

<v 2  >  „  -26  <V>2 

'vert  '  [ ln(Z/Zo ) ) 2 

where  Z0  ■  a  surface  roughness  scale  of  approximately  ore  to  five 
centimeters 

Z  *  the  height  of  the  laser  beam  above  ground. 

The  final  approximate  relationship  for  <6i>  is 

<62>  =  .28/[cos(4>)  In(Z/Z„) ] 2 . 


2 


The  result  is  that  <cos2(0)>  =1.0  and  <sin2(0)>  =  <02>  -  r  which 
are  then  used  to  derive  (2-9). 

The  major  restrictions  arising  from  the  above  approximations 
and  assumptions  are  that  n  be  very  small  and  that  angle  <J>  be  less 
than  plus  or  minus  45  degrees  which  restricts  angle  0  to  fairly 
small  values. 

Greenwood's  expression  for  the  spectrum  exhibits  the  sane 
-8/3  and  -2/3  power  dependencies  at  high  and  medium  frequencies 
as  did  the  spectra  derived  by  Clifford  and  Collins.  In  the  low 
frequency  region,  however,  only  Greenwood's  spectrum  has  a  flat¬ 
tening  trend  for  small  orientation  angles.  This  low  frequency 
flattening  trend  was  typical  of  experimental  data  measured  under 
conditions  of  strong  to  medium  turbulence  as  is  shown  by  our 
experimental  spectra. 
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Chapter  3 

DESCRIPTION  OF  EXPERIMENT 


3 . 1  Introduction 

Optical  phase  difference  experiments  at  10.6  micror  wave¬ 
lengths  were  conducted  from  February  through  May,  1972  at  the  Rome 
Air  Development  Center  Verona  test  site.  The  following  is  a  descrip¬ 
tion  of  the  simultaneous  measurements  that  were  made  of  optical  and 
meteorological  parameters. 

3 • 2  Optical  Measurements 

For  the  optical  measurements,  laser  beams  were  propagated 
through  the  atmosphere  for  a  distance  of  300  meters  over  a  path 
parallel  to  and  one  meter  above  the  ground.  The  beam  was  propagated 
between  the  transmitter  and  receiver  buildings  in  a  direction  36° 
south  of  west  over  uniform,  grassy,  level  terrain.  The  prevailing 
winds  were  from  the  northwest  at  an  angle  of  80°  to  the  propagating 
laser  beams.  However,  the  prevailing  winds  were  not  always  present 
on  days  that  experiments  were  conducted  as  will  be  shown  by  the  data. 

The  transmitter  was.  located  in  a  small,  weather-proofed  shelter 
and  mounted  on  a  stable  concrete  slab.  The  transmitter  consisted 
of  a  7-watt,  continuous-wave,  10.6-micron  wavelength  carbon  dioxide 
laser  and  a  50-milliwatt,  CL6328-micron  helium-neon  laser  whose 
beams  were  combined  and  transmitted  coaxially.  Before  leaving  the 
transmitter,  the  two  beams  were  demagnified  to  approximate  a  spherical 
wave  source.  The  helium-neon  beam  was  used  only  for  pointing  and 
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directing  the  carbon  dioxide  beam  into  the  receiver.  The  measure¬ 
ments  were  performed  upon  the  carbon  dioxide  laser  beam. 

3.3  Optical  Receiver 

The  optical  receiver,  shown  schematically  in  Fig.  3.1,  consisted 
of  a  16-inch  diameter  Boiler  and  Chivens  Cassegrain  telescope  and  an 
optical  propagation  effects  analyzer  built  by  the  Perkin-Elmer  Company. 

The  telescope  collected  the  energy  of  the  received  wavefront  and 
focused  it  at  the  primary  focus  of  the  optical  analyzer.  The  tele¬ 
scope  was  of  astronomical  quality  with  optics  of  one-tenth  of  a 
wavelength  smoothness  at  0.6  microns.  Consequently,  effects  from 
the  smoothness  of  the  telescope  surfaces  were  negligible  for  our 
experiments  at  10.6  microns.  The  Perkin-Elmer  optical  propagation 
analyzer  was  capable  of  performing  various  optical  light  beam  experi¬ 
ments.  For  this  experiment,  the  analyzer  ,/as  placed  in  a  twin- 
aperture-interferometer  configuration  to  measure  the  optical  phase 
difference  between  two  projected  sampling  aperture  points  on  the 
received  wavefront.  For  convenience  in  the  optical  analyzer  block 
diagram,  Fig.  3.1,  the  Cassegrain  receiving  telescope  and  the  two 
reflective  Dall-Kirkham  systems  are  illustrated  as  refractive  objectives. 
As  shown  in  Fig.  3.1,  the  primary  Dall-Kirkham  of  the  analyzer  accom¬ 
plished  an  eight  to  one  reduction  in  size  of  the  wavefront  received 
by  the  telescope.  The  two  sampling  apertures  of  the  interferometer 
were  positioned  in  the  collimated,  reduced  output  of  the  primary 
Dall-Kirkham.  The  secondary  Dall-Kirkham  converged  the  illumination 


of  the  two  apertures  and  produced  an  interference  pattern  at  its 
focal  point  which  was  on  the  surface  of  a  reflective  chopping  reticle. 


3 . 4  Phase  Measurement  Technique 

Appendix  A  contains  an  analysis  of  a  typical  twin  aperture 
interferometer  which  shows  that  the  intensity  distribution,  T(x,y,t), 
of  the  two-dimensional  interference  pattern  was 

cos  f K  x+4>  (t)  ]  }  (3-1) 


where  $(t)  is  the  time  varying  optical  phase  difference  between  the 


two  apertures 
N  -  4ttD2  (d/2)  4 
Kt  »  (TTd/Xf) 

K2  -  (2na/\f) 

a  is  the  aperture  separation 
d/2  is  the  aperture  radius 

D  is  the  amplitude  of  the  electric  field  at  each  aperture 
f  is  the  focal  length  of  the  convergent  optical  system. 

The  parameters  of  the  interferometer  were  such  that  K2>K1.  This  resulted 
in  an  interference  pattern  consisting  of  a  large  scale  circularly  sym¬ 
metrical  Airy  disc  pattern,  { [ 2 J x  (Kj  (x2+y2) ly^)  ] / [Kj  (x2+y2)  ^ ] }  which 
was  modulated  in  the  x-direction  by  the  higher  frequency  interference 
fringes  of  the  term  [1  +  cos(K2x+$(t) )  ] .  Equation  (3-1)  for  I(x,y,t) 
shows  that  the  position  of  the  fringes  in  the  x-direction  was  propor¬ 
tional  to,  4> (t ) ,  the  opticai  phase  difference  of  the  sampling  apertures. 
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In  the  optical  analyzer,  the  focal  point  of  the  interferometer 
and  hence  formation  of  the  interference  pattern  occurred  at  the 
surface  of  an  optical  reticle  which  was  rotating  at  a  constant 
angular  velocity  of  w  radians  per  second.  The  reticle,  illustrated 
in  Fig.  3.2,  was  a  quartz  disc  approximately  eight  inches  in  diameter. 
The  front  surface  of  the  reticle  consisted  of  a  periodic  pattern  of 
alternately  reflecting  and  absorbing  slits  of  spatial  frequency  K2 

and  of  radial  length  L.  As  shown  in  Fig.  3.2,  the  interference 

/ 

pattern  and  reticle  were  arranged  so  that  the  x-directicn  of  the 
interference  pattern  was  aligned  with  the  direction  tangential  to 
the  edge  of  the  reticle. 

Rotation  of  the  reticle  produced  a  modulation  of  the  inter¬ 
ference  pattern  intensity  which  was  reflected  from  the  reflective 
reticle  slits  to  a  mercury-cadrainn-tellur ide  detector.  This  modu¬ 
lation  was  equivalent  to  multiplication  of  the  Interference  pattern 
intensity,  J(x,y,t),  by  a  reticle  function,  M(x,y,t),  to  produce 
the  resultant  intensity 

I j (x,y,t)  -  I(x,y,t)M(x,y,t) .  (3-2) 

The  reticle  function  is  approximately  an  infinite  series  of  square 
waves  in  the  x-direction  of  spatial  frequency  K2  which  were  spatially 
phase  modulated  in  the  x-direction  :t  a  rate  corresponding  to  the 
tangential  velocity  of  the  reticle..  Thus,  for  slits  of  radial 
length  L  in  the  y-direction,  the  approximate  reticle  function 
becomes 

M(x,y,t)  -  [sQ(K2x  +  Wt)  [u(y+L/2)  -U(y-L/2)]]  (3-3) 
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where  SQ(K2x)  is  the  spatial  square  wave  function  of  frequency  K2 
R  =  radial  distance  on  reticle,  surface  to  the  center  of  the 
interference  pattern 

w  =  constant  angular  velocity  of  the  reticle  in  radians 
per  second 

U(y)  =  unit  step  function  in  the  y-direction 
W  =  K2Rw. 

Appendix  B  contains  an  analysis  of  the  error  that  arises  from 
approximating  the  reticle  function  by  an  infinite  series  of  square 
waves.  There  are  two  errors  which  occur  for  each  fringe  in  the 
interference  pattern.  The  first  error  is  in  the  location  or  absolute 
position  of  the  spatial  interval  corresponding  to  that  particular 
2tt  radian  fringe  period,  and  the  second  error  is  in  the  spatial 
length  of  each  fringe  period.  For  the  central  fringe  in  the 
pattern,  the  positional  and  interval  length  errors  ere  equal  and 
are  7.44  x  10~4  percent  of  one  interval.  This  is  much  smaller  than 
the  minimum  resolution  of  our  phase  measurement  which  was  one  part 
in  256.  Due  to  the  power  of  the  transmitter,  the  noise  equivalent 
power  level  of  the  infrared  detector,  and  the  input  noise  voltage 
level  of  the  detector  preamplifier;  the  power  level  of  the  seventh 
fringe  was  at  all  times  below  the  minimum  detectable  signal  level. 

The  errors  for  the  seventh  fringe  therefore  will  be  the  largest 
to  effect  the  resolution  of  the  experiment.  The  positional  error 
for  the  seventh  fringe  is  0.511  percent  of  one  interval  and  the 
interval  length  error  is  0.189  percent  of  one  interval.  Therefore, 
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in  the  worst  case  the  phase  resolution  is  limited  to  0.701  percent 
of  one  interval  or  O.OM*  radians. 

The  power  received  by  the  infrared  detector  is  the  integral 
over  the.  detector’s  field  of  view  (FOV)  of  the  product  of  the  modulated 
intensity  Ij(x,y,t)  and  one-half  the  admittance  of  free  space,  Y(.  The 
instantaneous  power  received  by  the  detector  was  therefore 

F(t)  »  -jY  //Ii  (x,y,t)  dydx. 

FOV 


Since  the  minimum  detectable  signal  level  limits  the  usable 
field  of  view  of  the  detector  in  the  interference  pattern  plane 
to  the  range  -B/2<x,y<+B/2 ,  the  expression  for  the  instantaneous 
power  becomes 


P(t) 


+B/2  +B/2 

±Y  N /  f 

2  0  -B/2  -B/2 

[SQ  (K2x+Wt) ] 


[  1+cos [K2x+  <f)(t)]  ] 


[U (y+L/2)  -  U (y-L/2) ]  dydx. 


(3-5) 


The  preceding  integral  can  be  simplified  by  the  following  assumptions 
which  apply  to  this  particular  case.  The  square  wave  function  can  be 
approximated  by  its  fundamental  term  cos (K2x+Wt) .  The  two-dimensional 
Airy  function,  ( [ 2J,  (Kx  (x2+y2)  ]/K,  (x2+y  *)  ^}2  ,  can  he  approximated 

by  a  constant.  A,  in  both  directions  since  the  effective  field  of  view 
of  the  detector  at  the  reticle  was  spatially  limited  to  a  small,  centrally 
located,  portion  of  the  Airy  function,  and  since  absolute  intensity  was 
not  of  importance. 

The  simplified  integral  for  the  Instantaneous  power  received  by 
the  detector  is 

+B/2  +B/2 

P(t)  -  jY  NA /  I  [  1+cos (Kzx+4>(t))]  [cos (K,x+Wt)] 

-B/2  -B/2  (3-6) 

[U (y+L/2 )-U (y-L/2 )]dydx 
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(3-7) 


which  after  the  y-integration  reduces  to 

j  +B/2 

P(t)  =  7Y0NAB  /  [l+cos(K2x+4>(t))][cos(K:2x+Wt)]  dx 

since  L/2>B/2. 

Assuming  the  field  of  view  in  the  x— direction  encompassed  an  integer 
number  of  cycles,  the  result  of  the  x- integration  is 

P(t)  =  jY#NAB  (2Tt/K2)  cos[<J>(t)-Wt]  (3-8) 

or 

P(t)  =  P  cos  [4>(t)-Wt ) .  (3-9) 

The  output  voltage  of  the  detectors  is  linearly  related  to  the 
input  power  by  the  detector  responsivity ,  R.  The  detector  output 
voltage  is  therefore 

V(t)  -  RP (t )  (3-10) 

or 

V(t)  *  RP  cos[<Kt)-Wt]  •  (3-11) 

These  equations  show  that  the  detector  input  power  and  output  voltage 
were  both  temporal  cosinusoidal  functions  phase  modulated  by  the 
optical  phase  difference  at  the  interferometer  aperture. 

A  phase  reference  signal  was  obtained  from  the  analyzer  by 
detecting  the  output  of  a  fixed  visible-wavelength  reference  light 
source,  the  intensity  of  which  was  also  modulated  by  the  chopping 
reticle.  The  resulting  reference  signal  voltage  was 

Vref(t>  “  1  cos[^ref-Wt].  (3-12) 

The  two  signals  from  the  optical  analyzer,  the  phase  measure¬ 
ment  signal  and  the  reference  signal,  were  then  processed  in  real 
time  by  a  hard-wired  digital  phase  difference  meter  which  measured 


31 


the  phase  difference  between  the  reference  channel  signal  and  the 
phase-modulated  measurement  channel  Gignal.  Fig.  3.3  is  a  siropji- 
fied  block  diagram  of  the  digital  phase  difference  meter.  The 
input  section  of  the  meter  performs  dual  channel  analog  signal 
conditioning  by  the  use  of  bandpass  filters  and  high  gain  compa¬ 
rators  with  hysteresis.  The  output  signals  of  the  bandpass  filters 

were  V  *  Pcos[4)(t)-Wt  ]  for  the  measurement  channel  and  Vrgf  * 
m 

Icos[<f  £ -Wt ]  for  the  reference  channel.  The  comparators  convert 
these  cosinusoidal  signals  into  fast  rise  time  square  waves  suitable 
fer  use  with  digital  circuits  for  the  following  digital  phase  differ¬ 
ence  measurement.  The  output  signals  from  the  comparators  are  the 
square  waves  Vm  -  S0[4>(t)-Wt]  and  =  SQ[$ref-WtJ  for  the  measure¬ 

ment  and  reference  channels,  respectively.  These  square  waves  are 
then  processed  by  two  negative-edge-triggered,  J-K  flip-frops  and 
an  F.xclusive-OR  gate  to  produce  a  pulse  whose  width  equals  the  time 
interval  corresponding  to  the  phase  difference  between  the  two 
channels.  This  pulse  then  controls  an  AND  gate,  permitting  a  parti¬ 
cular  number  of  clock  pulses  to  propagate  through  the  AND  gate  to 
the  phase  measurement  counters.  Thus,  the  number  of  caock  pul3es 
reaching  the  phase  measurement  counters  corresponds  to  the  Exclusive- 
OR  gate  pulse  width  which  is  equal  to  [ 4>(t)-4>ref  ] ,  the  difference 
between  the  time  varying  optical  phase  difference  and  a  constant 
arbitrary  phase  reference. 

At  this  point  the  phase  meter  performs  two  different  simul¬ 
taneous  phase  measurements;  one  of  the  intra-cycle,  within  2*  radians 
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instrumentation 


phase  difference  and  one  of  the  inter-cycle,  greater  than  2n  radians, 
phase  difference.  The  intra-cycle  phase  measurement  was  performed 
by  a  modulo-256  counter  and  an  8-bit  digital  to  analog  converter. 

The  resolution  of  the  intra-cycle  phase  measurement  was  eight  bits. 

The  inter-cycle  phase  measurement  was  performed  by  a  more 
complex  method  shown  schematically  in  Fig.  3.4.  The  clock  pulses 
representing  the  phase  difference  are  counted  by  a  modulo-256  counter 
whose  output  is  converted  from  digital  to  analog  form  by  an  8-bit 
converter.  The  differential  amplifier  takes  the  difference  between 
its  two  input  analog  voltages,  one  representing  the  current  phase 
difference  count  and  the  other  representing  the  previous  phase 
difference  count  which  was  stored  by  the  sample-and-hold  circuit. 

The  difference  between  the  present  and  previous  phase  counts  is 
compared  to  a  threshold  level  by  the  analog  comparator.  When  the 
phase  difference  exceeds  2iT  radians  the  input  to  the  modulo-256 
counter  exceeds  256,  the  2^  level,  and  the  counter  continues 
counting  the  phase  difference  starting  from  zero  again.  Due  to 
the  rapid  sample  rate  at  which  the  measurement  is  performed,  in 
one  cycle  for  the  reference  channel  the  measurement  channel  phase 
cannot  change  more  than  several  degrees.  This  results  in  a  large 
difference  between  the  present  and  previous  phase  measurements  when 
a  slip  of  one  cycle  occurs.  This  difference  exceeds  the  threshold 
level  of  the  comparator  and  causes  its  output  to  change  state. 

The  comparator  output  is  then  differentiated  by  a  high  pass  filter 
to  form  a  fairly  sharp  pulse  which  in  turn  causes  the  output  up-down 


counter  to  count  up  or  down  by  one,  the  direction  of  the  count 
dependent  upon  the  direction  of  the  phase  slip.  This  counter  is 
not  reset  at  each  sample  and  hence  represents  the  long  term  inter¬ 
cycle  phase  difference.  The  output  counter  is  also  connected  to  a 
digital-fo-analog  converter  which  has  a  resolution  of  six  bits.  When 
propagation  conditions  caused  the  dynamic  range  of  the  phase  measure¬ 
ment  to  increase  beyond  2tt  radians,  it  was  then  possible  to  combine 
these  two  separate  measurements  to  yield  the  required  dynamic  range. 
3. 5  Meteorological  Measurements 

Simultaneous  to  the  performance  of  the  optical  experiments, 
measurements  were  made  of  gross  meteorological  parameters  and  of  the 
microstructure  of  the  temperature  field. 

The  wind  direction  and  speed  are  important  in  determining 
scaling  parameters  for  analyzing  optical  data.  Wind  direction  was 
measured  with  a  Packard-Bell  low-inertia  wind  vane  system  which  had 
a  -3dB  point  at  5Hz.  Wind  direction  measurements  were  made  with 
Packard  Bell  six-cup  anemometers  which  had  a  -3dB  point  also  at  5Hz. 

Micro  temperature  was  measured  on  a  single  probe  constructed 

(9) 

at  RADC  according  to  the  design  of  Ochs  at  NOAA.  The  RADC  sen¬ 
sors  are  platinum  wire  2.5pm  in  diameter  and  about  2  mm  in  length. 

The  3db  point  of  such  a  probe  is  about  500  Hz  in  5m/sec  winds,  and 
that  is  certainly  adequate  for  expected  fluctuation  rates.  The 
sensor  is  connected  as  one  leg  of  a  Wheatstone  Bridge.  Tempera¬ 
ture  changes  sensed  by  the  probe  cause  a  proportional  imbalance 
voltage  at  the  bridge.  The  imbalance  voltage  is  amplified  by  a  gain 
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of  1.1  x  10  and  then  recorded  on  analog  tape  for  later  processing.  As 
the  orobe  draws  only  about  450pa  current,  it  operates  as  a  cold-wire 
sensor  for  wind  speeds  exceeding  0.5  m/sec.  A  much  more  extensive 
description  of  the  micror hermal  system  is  given  in  a  report  by  Greenwood 
and  Tarazano.^^ 

3 . 6  Data  Recording  and  Reproduction 

Wien  the  experiments  were  conducted,  the  two  channels  of  optical 
phase  difference  data  from  the  digital  phase  meter  and  all  the  meteo¬ 
rological  data  were  recorded  on  an  Ampex  analog  instrumentation  tape 
recorder.  The  frequency-modulation  mode  of  recording  was  used  at 
tape  speeds  of  eitner  7.5  or  15  inches  per  second  which  yielded 
recorder  bandwidths  of  0  Hz  to  2500  Hz  and  0  Hz  to  5000  Hz,  respectively. 

The  data  was  later  replayed  on  the  tape  recorder  tor  digiti¬ 
zation  by  a  Digital  Equipment  Corporation  PDP-8  computer  with  a  12-bit 
analog-to-digital  converter.  Prior  to  digitization,  the  signals  were 
filtered  to  prevent  aliasing  with  four-pole  Butterworth  low-pass 
filters  set  for  the  proper  antialiasing  at  the  Nyquist  frequency. 


Chapter  4 


ANALYSIS  OF  EXPERIMENTAL  DATA 


4.1  Introduction 

The  optical  and  meteorological  data  collected  during  the 
experiment  was  initially  recorded  in  analog  form.  The  data 
recordings  were  replayed  at  a  later  time  for  statistical  and  power 
spectral  density  analysis. 

4.2  Optical  Data  Analysis 

The  optical  phase  difference  data  was  processed  to  yield 
estimates  of  the  average  temporal  power  spectral  density  of  the 
phase  difference  fluctuations.  Most  of  the  optical  power  spectral 
density  analysis  and  processing  was  performed  by  a  fast  Fourier 
transform  technique  on  a  digital  computer.  However,  some  of  the 
high  frequency  power  spectral  density  information  was  obtained 
by  processing  the  data  with  a  Ubiquitous  Spectrum  Analyzer  and 
Averager  system.  This  chapter  wfll  describe  the  analysis  pro¬ 
cedures  and  the  associated  confidence  limits  of  the  resulting  data 
for  both  types  of  power  spectral  density  processing. 

For  the  fast  Fourier  transform  processing,  the  analog  data 
was  converted  to  12-bit  digital  words  and  then  stored  on  digital 
magnetic  tapes  by  a  Digital  Equipment  Corporation  PDP-8  computer 
equipped  with  a  12-bit  analog-to-d-'gital  converter.  These  digital 
data  tapes  were  then  processed  by  a  Fortran  TV  fast  Fourier  transform 


The  Honeywell  computer  used 


algorithm  on  a  Honeywell-645  computer, 
a  16-bit  floating-point  arithmetic  format  which  had  27  binary  bit 
significance  in  the  mantissa  of  the  floating-point  woids.  There 
are  three  major  factors  affecting  the  accuracy  of  the  processed 
P°v'er  spectral  density  daca.  These  accuracy  determining  factors 
are  errors  that  arise  from  the  digitization  and  computer  processing 
of  the  original  analog  data. 

The  three  errors  in  the  fast  Fourier  transform  algorithm  pro¬ 
cessing  all  arise  from  computer  and  anolog-to-digital  converter 
roundtng  of  numerical  values  of  infinite  precision  to  values  with 
finite  precision.  The  error  due  to  rounding  by  the  12-bit  analog¬ 
ic-digital  conversion  of  the  initial  data,  x(nT),  may  be  represented 
by  an  additive  noise,  e(nT),  which  has  a  uniform  probability  distri¬ 
bution,  zero  mean,  variance  of  one-twelth  of  one  bit,  and  which  is 
uncorrelated  between  samples.  The  input  to  the  algorithm  is  then 
x(nT)  *  x0(nl)  -f  e(nT)  where  x0(nT)  is  the  noise  free  input  rounded 
to  12  bits.  The  evaluation  of  the  propagation  of  this  error  through 
the  fast  Fourier  transform  algorithm  by  analytical  means  is  a  diffi¬ 
cult  task.  An  experimental  test  of  the  effects  of  this  error  on 
the  dynamic  range  of  the  output  power  spectra  can  be  conducted  using 
a  "white  noise"  generator  and  filters  of  known  characteristics.  This 
test  was  not  conducted.  However,  the  experimental  data  which  was 
random,  Gaussian  distributed,  and  not  white  but  smoothly  filtered, 
typically  exhibited  a  linear  dynamic  range  of  63  dB. 

The  remaining  two  computer  generated  errors  are  due  to  rounding 
all  multiplication  products  and  multiplication  coefficients  to  produce 
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27-bit  mantissas.  For  a  typical  radix-two  fast  Fourier  transform 
algorithm  performed  with  floating  point  arithmetic,  the  variance 
of  the  error  in  the  transformed  data  which  arises  from  the  rounding 
of  the  multiplication  products  has  been  shown  to  be 

0  E2,RO  *  2/3  v  N  2  2b  Ox2  (4-1) 

where  N  a  4096,  the  number  of  points  in  the  input  sequence 
v  -  log2 (N) 

b  *  27,  the  number  of  bits  in  the  mantissa  of  the  floating 
point  word 

0x  «  the  variance  of  the  input  sequence. 

Since  for  the  experimental  data  the  maximum  input  signal  variance 
was  3.2  x  10+J  bits  squared,  the  maximum  root-mean-squat e  output 
error  due  to  rounding  of  multiplication  products  was  approximately 
0.008  bits. 


The  variance  of  the  error  in  the  transformed  data  due  to 


rounding  of  multiplication  coefficients  has  been  shown 


^  to  be 
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which  for  the  case  of  N  ■  4096  reduces  to 


E,C 


**  ta 


E,R0  . 


(4-3) 


The  root-mean-square  error  in  the  transformed  sequence  due  to 
coefficient  rounding  is  approximately  0.002  bits.  Therefore, 
considering  both  multiplication  product  rounding  and  coefficient 
rounding,  the  root-mean-square  output  error  was  less  than  one  bit. 

In  summary,  for  the  fast  Fourier  transform  processing,  the 
analog-to-digital  converter  round-off  error  was  experimentally 
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determined  to  limit  the  dynamic  range  of  the  processed  spectra  to 
60  dB,  or  approximately  ten  bits.  The  error  generated  by  the  rounding 
of  both  the  multiplication  products  and  the  multiplication  coefficients 
in  the  algorithm  was  shown  to  introduce  a  root-mean-square  error  of 
less  than  one  bit  in  the  27-bit  mantissas  of  the  transformed  data. 
Therefore,  the  analyzed  spectra  have  a  dynamic  range  limited  to 
approximately  60  dB  or  ten  bits  by  propagation  of  the  analog-to- 
digital  converter  round-off  error  through  the  algorithm. 

Some  of  the  experimental  power  spectral  density  data  showed 
a  large  spread  between  individual  spectral  curves  at  low  frequencies. 
As  discussed  in  chapter  5  this  spectral  spread  is  due  to  the  large, 
random  wind  direction  fluctuations.  To  conclude  that  the  spread 


of  the  low  frequency  spectra  was  due  to  the  wind  direction  fluctu 
ations,  it  was  helpful  to  derive  appropriate  confidence  or  expec¬ 


tation  intervals  for  the  spectra  based  on  the  assumptions  of  station- 
arity  and  a  chi-square  distribution  for  averaged  spectral  values. 

The  confidence  limits  are  due  to  the  fact  that  the  resulting  pro¬ 
cessed  power  spectral  density  values,  P*^0)»  arc  ranc*o:D  variable 
estimates  of  the  true  power  spectral  density,  W(fQ),  of  the  initial 
process.  Each  power  spectral  density  estimate,  P^(fQ),  is  actually 
a  random  variable  defined  as  the  sample  mean  of  a  set  of  discrete. 


al  density  value  estimate. 


p2(f  )  is  the  sample  mean  power  spectr 
B“  making  suitable  assumptions  It  Is  possible  to  develop  an 
jpproxlmate ,  ganaral  ctltetlon  lot  the  confidence  limit*  of  the  estimate 
ahlch  docs  not  require  knowledge  of  the  statistics  of  the  original  pro 
ress.  The.  original  process  Is  assumed  to  be  stationary,  Gaussian  dis 
tributed,  and  to  possess  a  relatively  smooth  power  speettum.  Hence, 
the  sample  mean  power  spectral  'density  estimates  will  be  chi-square 

distributed  random  variables. 

The  probability  statement  that  detincs  the  interval  about 

P*(f  ),  v,  the  sample  near,  estimate,  Ln  which  we  will  find 

or  theoretical  power  spectral  density  value,  W(fQ),  with  c  proba 

of  0.8  is 

Prob  t'Y/cl0)<(K(fo))<(^o0)l  "  °'8  ^ 

,  <-  f-Mvi  nrobabilitv  statement  bv  showing 

Fig.  5.2  illustrates  the  use  of  this  niobaoii  . 

i  mrve  and  its  associated  80 

one  experimental  power  spectral  density  curve  ana 

percent,  rontidence  intervals.  The  sire  of  the  confidence  interval 
decreases  for  higher  frequencies  because  a  larger  number  of  individual 
po„er  spectral  values  were  used  ro  for.  the  sample  mean  estimate. 

The  probability  statement  defining  the  interval  about  «<»„>.  the 

i  .  4-  which  we  expect  to  find  80  percent 

true  power  spectral  density  value,  ...  which 

of  the  power  spectral  density  estimates,  V^t Q) ,  *■* 

Prob  t(C90W(f0))<Y<<c»oH(fo))1  "  °'8, 

^  i  mn-v  R«-"tement  is  illustrated  In  Figures  '>.3 

Thf.  use  of  this  probability  st-tenu-m. 

.  about  the  theoretical 

and  5,6  where  the  80  pec  ent  expectation  intervals 

power  spectral  curve  are  shown. 


(4-6) 
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For  the  10  percent  and  90  percent  confidence  limits  with  which 

(8) 

we  are  concerned,  che  quantity  C  is  tabulated  for  n  <  30.  The 

P 

equations  given  in  this  reference  show  that  suitable  approximations 

for  C  when  n  >  30  are 
P 

C10  =  1.0  +  (1. 81/n/n)  (4-7) 

C90  “  1.0  -  (1.8lA/n)  (4-8) 

Tt  should  be  noted  that  the  preceding  assumptions  of  statiori- 
arity  and  a  chi-square  distribution  of  power  spectral  density  esti¬ 
mates  cause  these  general  case  confidence  Limits  to  be  approximations 
whose  applicability  depends  upon  the  degree  of  adherence  to  the  above 
restrictions.  Hence,  for  some  of  the  data  which  is  widely  distributed 
without  a  chi-square  distribution  of  power  spectral  values,  the  confi¬ 
dence  limits  indicated  are  not  representative  of  the  true  sample  mean 
and  variance 

The  fast  Fourier  transform  and  power  spectral  density  computer 
program  calculates  2,048  discrete  power  spectral  values.  The  discrete 
spectral  values  occur  at  intervals  of  1/4,096  of  the  sample  rate  over 
the  range  of  zero  hertz  to  one-haif  of  the  sample  rate.  The  computer 
program  then  condenses  these  2,048  spectral  values  into  a  set  of  32 
new  values  which  are  equally  spaced  along  a  logarithmic  abcissa  or 
frequency  axis.  These  new  values  are  formed  by  averaging  an  appro¬ 
priate  number  the  initial  2 , C48  values  in  the  regions  centered 
at  e3ch  of  the  32  values.  Table  4.1  lists  the  center  frequencies 
corresponding  to  each  of  the  32  discrete  intc’-vals  for  tue  sample 
frequency  of  ten  herti  which  was  used  for  the  low  frequency  data 
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Table  4.1 


Data  Pertaining  to  the  Low-frequency  Power 
Spectral  Density  Curves 


Center 

Frequency 


0.0024 
0. J049 
0.0073 
0.0098 
0.0122 
0.0146 
0.0171 
0.0220 
0.0293 
0.0366 
0.0464 
0.0586 
0.0732 
0.0903 
0.1123 
0.1416 
0.1782 
0.2222 
0.2783 
0.3491 
0.4370 
0.5469 
0.6812 
0.8521 
1.0596 
1.3159 
1.6455 
2.0605 
2.5732 
3.2202 
4.0259 
4.7339 


Number  of 
Points  Averaged 


1 

1 

1 

1 

1 

1 

1 

3 

3 

3 

5 

5 

7 

7 

11 

13 

17 

19 

27 

31 

41 

49 

61 

>9 

°1 

119 

151 

189 

231 

399 

361 

218 


80%  Confidence  Limit 


Factors 


Cso 

Ci  o 

0.0158 

2.706 

0.0158 

2.706 

0.0158 

2.706 

0.0158 

2.706 

0.0158 

2.706 

0.0158 

2.706 

0.0158 

2.706 

0.1967 

2.084 

0.1967 

2.084 

0.1967 

2.084 

0.322 

1.847 

0.322 

1.847 

0.405 

1.717 

0.405 

1.717 

0.507 

1.57 

0.542 

1.52 

0.593 

1.46 

0.613 

1.43 

0.671 

1.36 

0.675 

1.33 

0.717 

1.28 

0.741 

1.26 

0.768 

1.23 

0.796 

1.20 

0.810 

1.19 

0.834 

1.17 

0.853 

1.15 

0.868 

1.13 

0.881 

1.19 

0.909 

1.09 

0.905 

1.10 

0.877 

1.12 
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The  table  also  lists  the  number  of  points  averaged  over  the  32  inter¬ 
vals  and  the  resulting  10  percent  and  90  percent  confidence  limits 
for  the  resulting  32  mean  spectral  values.  The  tabulated  confidence 
limit  values  are  the  factors  by  which  the  experimental,  mean  spectral 
values  would  be  divided  to  yield  the  range  in  which  we  have  80  percent 
confidence  of  finding  the  true  or  theoretical  power  spectral  density 
values . 

To  reduce  the  large  confidence  limits  which  exist  at  the  lowest 
frequencies,  individual  power  spectral  curves  were  averaged  together. 
These  final  mean  power  spectral  curves  for  each  day  of  experimental 
data  are.  presented  in  chapter  5  and  error  bars  for  the  80  percent 
confidence  limits  are  also  shown. 

The  optical  high  frequency  data  was  analyzed  with  a  Federal 
Scientific  Corp.  Ubiquitous  Spectrum  Analyzer  (Model  UA-10A)  and 
Averager  (Model  1010).  The  Ubiquitous  Analyzer  performs  a  direct 
power  spectral  measurement  by  means  of  a  narrow  bandpass  filter  and 
power  detector  technique.  However,  the  Ubiquitous  Analyzer  had  a 
linear  dynamic  range  of  only  A0  dB.  Therefore,  the  data  had  to  be 
divided  into  segments  whose  spectrum  did  not  exceed  this  AO  dB  range. 
For  this  reason  only  the  A0  dB  high-frequency  portion  of  the  spectrum 
was  processed  with  the  Ubiquitous  Analyzer.  It  was  then  possible  to 
process  the  entire  lower  portion  of  the  data's  spectrum  with  the 
computer  fast  Fourier  transform  techniques  which  had  a  60  dB  dynamic 
range.  The  high-frequency  spectral  measurements  were  performed  by 
averaging  512  individual  spectra  to  produce  the  final  mean  spectrum. 
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The  resulting  mean  spectra  had,  therefore,  512  degrees  of  freedom  and 
the  corresponding  upper  and  lower  confidence  limits  are  1.08  and  0.9? 
times  the  spectral  values.  The  corresponding  error  bars  are  relatively 
small  and  have  not  been  shown  on  the  spectral  curves. 

4.3  Meteorological  Data  Anal  vs Is 

The  meteorological  data  processing  consisted  of  computing  the  sample 
mean  and  sample  variance  of  the  recorded  wind  speed  and  wind  direction  data, 
and  computing  the  power  spectra  of  the  microtemperature  fluctuations. 

The  processing  of  the  wind  speed  data  was  straight  forward. 

The  data  was  replayed  from  the  analog  tape  recordings  and  digitized 
by  the  PDP-8  computer  at  a  sample  rate  of  approximately  ten  hertz. 

The  wind  speed  systems  have  a  frequency  response  of  approximately 
5.0  Hz  which  was  used  as  the  anti-aliasing  filter.  Real  time  calcu¬ 
lation  of  the  sample  mean,  n,  and  sample  variance,  o2,  was  performed 
with  the  computer  according  to  the  formulas 

n  -  ^mlXt  (4-9) 

U2.ifi(X2)  (4-10) 

Ni-l  i 

-  u2-  n2  (4-11) 

where  M  is  the  second  moment  of  the  sampled  data. 

The  wind  direction  data  processing  was  performed  in  the  same 
manner  as  was  the  wind  speed  processing  except  that  a  different  com¬ 
puter  program  was  required  to  perform  the  calculations.  The  wind 
direction  transducer  represented  wind  directions  of  0°  to  360°  by  an 


analog  output  voltage  of  zero  to  five  volts.  Therefore,  when  the  wind 
direction  fluctuated  about  the  0°,  360°  crossover  point,  the  output 
voltage  fluctuated  from  zero  to  five  volts.  If  the  previous  formulas 
for  sample  mean  and  variance  were  used  for  this  type  of  data,  incorrect 
values  for  n  and  c  would  result.  Due  to  this  anomaly,  a  special  pro¬ 
gram  was  required  to  calculate  the  mean  wind  direction. 

The  mean  wind  direction  program  sectioned  the  digitized  wind 
direction  data  into  the  four  quadrants  of  the  compass  circle,  counted 
the  number  of  samples  per  quadrant  and  calculated  the  mean  wind  direction 
per  quadrant .  To  better  characterize  the  magnitude  of  wind  direction 
fluctuation,  values  for  the  mean  angle  per  quadrant  and  the  number  of 
samples  per  quadrant  were  calculated  and  displayed  nine  times  per 
data  record.  At  the  end  of  the  data  record  the  mean  angle  for  the  entire 
record  was  computed  by  transforming  the  mean  angles  for  each  quadrant  to 
their  rectangular  components  along  the  quadrant  axis,  weighting  these 
components  by  the  number  of  samples  Per  quadrant,  and  computing  the  mean 
angle  trom  their  sum.  This  also  permitted  characterization  of  the  frequency 
of  the  wind  direction  fluctuations  in  a  crude  fashion.  A  typical  computer 

printout  for  a  day  which  had  large  wind  direction  fluctuations  is  shown 
in  Table  5.2. 


The  power  spectral  analysis  o'  the  microthermal  signals  orllUed 
the  fast  rentier  transform  approach  described  In  Section  4.2.  The  major 
difference  Is  that  the  Ubiquitous  Spectrum  Analyser  was  used  in  the  high 
frequencies  there,  but  the  digital  computer  is  used  for  the  high  frequency 
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results  in  the  microthermal  analysis.  The  number  of  points  averaged 
together  for  the  high  frequencies  is  64  times  the  tabulated  values  in 
Table  4.1.  The  center  frequencies  are  also  shifted  by  a  factor  of  40. 
The  resultant  confidence  limits  may  be  easily  calculated  via  equations 
(4-7)  and  (4-8). 
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Chapter  5 
EXPERIMENTAL  DATA 

5*1  Introduction 

EXP'rl”"t*1  "ea— "»  •*  the  optical  phase  difference  fM 
'O’6'"1"0"  “6th  aad  ‘he  corresponding  meteorological 

performed  ^  dflys  during  the  monthg  ^ 

AP;U  “*  ^  *  1972‘  -  ~  „f  experimental  prohlems  ’ 

en  unfavorable  weather  oonditiona  there  remained  only  five  days  of 

acceptable  data.  The  experimental  data  of  two  typical  days  was  com- 
Pletely  analyzed  and  will  be  presented. 

The  initial  an.ly.is  of  the  optical  data  .bowed  that  the  low- 
frequency  region  of  the  phase  difference  power  spectrum  behaved 

differently  on  day.  of  erratic  wind  direction  condition,  than  on  days 
Of  steady  wind  direction.  Therefore  d.t.  t 

,  ,  ’  data  from  d*ys  °f  both  conditions 

is  Included. 


5,2  Experimental  Data 

Optical  and  meteorological  measurements  were  performed  on  the 

°f  May>  1,72  b”“""  d'dloch  in  the  morning  .„d  two 

°  Clock  In  the  afternoon.  This  was  a  clear  dav  with 

•Lear  day  with  a  mean  temperature 

of  76°F. 

Table  5.1  lists  the  means  and  variances  of  the  horizontal  wind 
speed  for  the  four  fourteen^inute  intervals  over  which  data  was 
recorded.  The  mean  wind  speed  fluctuated  from  1.11  meters  per  second 

10  P"  8CC°nd  “Uh  “  average  of  1.6  meters  per 
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second.  Thus  on  this  day  the.  wind  speed  was  very  low  and  had  signi¬ 
ficant  amount  of  variance.  Table  5.2  of  wind  direction  data  for  a 
typical  measurement  interval  shows  that  the  wind  direction  on  this  day 
had  large  low  frequency  fluctuations  greater  than  270°  in  magnitude. 

As  shown,  at  times  the  wind  direction  varied  slowly  back  and  forth 
across  the  laser  beam  propagation  path. 

Fig.  5.1  shews  the  complete  set  of  phase  difference  power 
spectra  measured  on  this  day.  There  are  eight  low-frequency  spectra 
each  measured  over  a  seven  minute  time  period  and  four  corresponding 
high-frequency  spectra  each  measured  over  two  of  the  seven  minute 
periods.  At  the  lower  frequencies  these  experimental  spectra  exhibit 
a  rising  trend  and  are  spread  over  four  decades  vertically.  This 
four  decade  vertical  spread  is  much  greater  than  the  spread  permitted 
by  the  80  percent  confidence  limits  of  these  spectra.  Hg.  5.2 
illustrates  this  by  showing  one  spectrum  from  Fig.  5.1  and  its  80 
percent  confidence  limits. 

Fig.  5  J  contains  a  graph  of  the  mean  experimental  spectrum, 
averaged  ever  all  the  individual  spectra  of  Fig.  5.1.  Fig.  5.3 
ait»o  contains  a  graph  or  Greenwood’s  theoretical  spectrum  plotted 
fer  an  outer  scale  size  of  two  meters,  a  sectnd  moment  of  the  orien¬ 
tation  angle  of  0.12  radians  squared,  and  a  mean  wind  speed  perpen¬ 
dicular  to  the  beam  of  0.6  meters  per  second.  A  mean  wind  speed  of 
0.6  metera  per  second  was  used  because  it  gave  the  best  fit  between 
the  theoretical  and  experimental  curves.  It  was  necessary  to  use 
this  method  because  the  large  variations  in  wind  direction  prohibited 
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Table  5.2 


Mean  Wind  Direction  Data  for  May  18,  1972; 
Data  Record  Number  A 


0°  to  90° 

90°  to  1808 

180°  to  270° 

270°  to 

Degrees  - 

AO. 6 

163.9 

26A.0 

286.0 

if  Samples  ■ 

1.0 

1.0 

1A3.0 

855.0 

Degrees  «■ 

9.9 

0.0 

268. A 

305.9 

if  Samples  ■ 

105 

0.0 

15.0 

88.0 

Degrees  - 

3.8 

0.0 

0.0 

333.7 

if  Samples  ■ 

282 

0.0 

0.0 

718.0 

Degrees  - 

72.6 

11A.1 

26A.0 

289. A 

if  Samples  ■ 

1.0 

2.0 

206.0 

791.0 

Degrees  - 

58.3 

133.1 

262.1 

27A.6 

if  Samples  - 

3.0 

A  .0 

A51.0 

5A2.0 

Degrees  - 

35.5 

119.  A 

261.5 

273.8 

if  Samples  ■ 

3.0 

7.0 

266.0 

72A.0 

Degrees  - 

0.0 

0.0 

250.9 

275.3 

if  Samples  ■ 

0.0 

0.0 

892.0 

108.0 

Degrees  » 

158.2 

116.  A 

2AA.A 

276. A 

if  Samples  ■ 

3.0 

2.0 

850.0 

1A5.0 

Degrees  ■ 

0.0 

0.0 

260.5 

292.1 

if  Samples  ■ 

0.0 

0.0 

A3.0 

9S7.0 
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FREQUENCY  (HERTZ) 

g.  5»1»  The  combined  set  of  high  and  low  frequency  experiment 
fference  spectra  measured  on  May  1 8,  1972. 


PRE^USNCY  (HERT2) 

Pig.  5.2.  An  individual  phase  difference  power  spectnun  mesaured  on  Key  IP,  1972. 
The  vertical  bars  indicate  the  60  percent  confidence  intervals  about  each  power 
spectral  value  (x  x  x). 
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2  m/s 


a  meaningful  determination  of  the  effective  perpendicular  wind  velocity 
along  the  entire  beam  from  a  single  point  wind  speed  and  wind  direction 
measurement . 

The  theoretical  and  mean  experimental  curves  fit  closely  together 
in  the  high-frequency  region  where  f>Y6Vj_/27Tp  or  f>3.2vL>  In  the 
medium-frequency  region,  defined  as  AvL/2irL0<f<  vJL/2irp  or  .08^  < 
f<3.2V|_,  the  two  curves  begin  to  separate  vertically  by  less  than  a 
decade.  This  separation  of  the  two  curves  at  medium  frequencies  is 
due  to  both  poorer  confidence  limits  and  the  same  condition  :  which 
caused  even  larger  separations  in  the  low  frequency  region.  In  the 
low-frequency  region  f<Av^ /2ttLo  ,  the  mean  experimental  curve  rises 
four  decades  above  the  flattening  level  of  the  theoretical  spectrum. 

The  micro temperature  power  spectra  for  May  18,  1973,  Figures  5.3a 
and  5.3b,  demonstrate  a  behavior  in  the  low-frequencies  similar  to  the 
phase-difference  spectrum.  The  data  in  figure  5.3b,  indicated  by  0's  are 
the  averages  of  the  individual  curves  in  Figure  5.3a.  The  theoretical 
curve  in  Figure  5.3b  is  based  on  von  Karman's  spatial  spectrum.  The 

excess  power  in  the  low-frequencies  of  the  microtemperature  spectrum 
insures  us  that  there  was  no  basic  error  in  the  measurement  of  phase- 
difference. 

The  separation  of  the  theoretical  and  experimental  spectra  at  low 
frequencies  is  too  great  to  be  entirely  accounted  for  by  the  existing 
confidence  limits.  Instead,  this  separation  must  be  attributed  to  the 
large  fluctuations  in  wind  direction  and  to  the  low  and  intermittent  wind 
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speed  conditions  which  existed  during  the  measurement  of  these  spectra. 
For  these  conditions  of  wind  speed  and  wind  direction,  the  large-scale 
•spatial  anisotropy  and  non-homogeneity  of  the  turbulence  and  refractive 
index  fields  contribute  a  significant  amount  of  variance  to  the  level  of 
the  spectra  at  low  frequencies. 

Consideration  of  the  local  terrain  as  sketched  in  Figure  5.4  shows 
the  existence  of  large  wooded  areas  along  the  north  side  of  the  propaga¬ 
tion  path  ana  immediately  behind  the  transmitter.  The  development  of  the 
aforementioned  anisotropic  and  non-homogeneous  turbulence  is  easily 
accounted  for  since  at  times  the  wind  was  from  the  wooded  areas  and  at 
other  times  from  the  clear  fields  south  of  the  path.  The  major  difference 
in  turbulence  generated  in  wooded  areas  as  opposed  to  turbulence  from 
unobstructed  level  terrain  is  at  low  spatial  frequencies  or  large  scale 
“lzc8‘  Therefore>  during  this  experiment  trie  fluctuating  wind  direction 
end  lew  wind  speed  conditions  produced  fluctuations  in  the  spatial  distri¬ 
bution  of  the  larger  scales  of  turbulence  which  were  optically  measured  as 
low  temporal  frequency  fluctuations  of  the  optical  phase  difference. 

The  second  set  of  experimental  measurements  were  performed  on  the 
25th  of  May,  1972  between  eleven  o’clock  in  the  morning  and  one  o'clock  in 

the  afternoon.  This  was  another  clear  day  for  which  the  mean  temperature 
was  69°F. 

The  horizontal  wind  speed  for  this  day  was  moderate  and  fairly 
steady.  Over  the  entire  data  recording  period,  the  mean  wind  speed  was 
5.5  meters  per  second.  Table  5.3  lists  the  mean  values  of  the  wind  speed 
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Table  5.3 


^ 


Wind  Speed  and  Direction  Data  for  May  25,  1972 


Data 

Record  Mean  Wind  Speed 
Number  (Metera/Second) 

1  6.14 

2  5.93 

3  5.20 

4  5.22 

5  5.32 


Mean  Wind  Direction 
Angle  (Degrees  from  North) 

82.7 
71.0 

70.8 

71.8 
77.5 


Mean  Horizontal 
Wind  Speed  Component 
Perpendicular  to 
Laser  Beam 
(Meters /Second) 

2.96 

1.65 

1.41 

1.73 

2.10 
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Table  5. A 


Mean  Wind  Direction  Data  for  May  25,  1972- 
Data  Decord  Number  A 


0®  to  90° 

90°  to  180° 

180°  to  270° 

270°  to  360 

Degrees 

ts 

70.5 

90.3 

0.0 

0.0 

it  Samples 

8A9.0 

151.0 

0.0 

0.0 

Degrees 

75.98 

0.0 

0.0 

0.0 

if  Samples 

ms 

1000.0 

0.0 

0.0 

0.0 

Degrees 

- 

67. 1A 

0.0 

0.0 

0.0 

if  Samples 

IS 

1000.0 

0.0 

0.0 

0.0 

Degrees 

m 

79.63 

0.0 

0.0 

0.0 

if  Samples 

MS 

1000.0 

0.0 

0.0 

0.0 

Degrees 

B 

78. A3 

0.0 

0.0 

0.0 

it  Samples 

B 

1000.0 

0.0 

0.0 

0.0 

Degrees 

- 

1A0.0 

0.0 

0.0 

0.0 

it  Samples 

m 

1000.0 

0.0 

0.0 

0.0 

Degrees 

= 

62.96 

0.0 

0.0 

0.0 

it  Samples 

t  * 

1000.0 

0.0 

0.0 

0.0 

Degrees 

- 

69.2 

90.1 

0.0 

0.0 

it  Samples 

m 

937.0 

63.0 

0.0 

0.0 

Degrees 

m 

6A.9A 

0.0 

0.0 

0.0 

it  Samples 

MS 

1000.0 

0.0 

0.0 

0.0 

FREQUENCY  (HERTZ) 

Pig.  5.5.  The  combined  set  of  high  and  low  freouenoy  experimental  phase 
difference  spectra  measured  on  May  25»  1972. 
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for  the  five  fourteen-minute  individual  data  records.  The  mean  wind 
direction  over  the  entire  data  recording  period  was  74.76°  from  the  north 
or  112.5°  from  perpendicular  to  the  path.  Table  5.4  which  contains  the 
wind  direction  data  for  a  typical  individual  data  record  shows  that  the 
wind  direction  variance  was  much  less  than  for  the  previous  data  of 
March  18th. 

Figure  5.5  shows  the  set  of  experimental  spectra  measured  on  May 
25,  1972.  The  ten  low-frequency  spectra  were  measured  over  a  seven-minute 
time  period  and  the  five  corresponding  high-frequency  spectra  were  each 
measured  over  two  of  the  seven-minute  periods.  As  clearly  illustrated  by 
the  graph,  seven  out  of  ten  of  the  spectra  exhibited  the  theoretically 
predicted  flattening  tendency  at  low  frequencies. 

Figure  5.6  shows  a  mean  experimental  spectrum  averaged  over  the 
seven  curves  which  had  a  flattening  trend  at  low  frequencies.  This  figure 
also  contains  a  graph  of  Greenwood’s  theoretical  spectrum  with  the  appro¬ 
priate  80  percent  confidence  limits.  This  theoretical  spectrum  is  for  an 
outer  scale  of  two  meters,  an  orientation  angle  variance  of  .12  radians 
squared,  and  for  a  mean  wind  speed  perpendicular  to  the  beam  of  1.0  meters 
per  second.  This  value  of  mean  perpendicular  wind  speed  is  one-half  the 
measured  value  of  1.97  meters  per  second,  but  it  is  the  perpendicular  wind 
speed  which  resulted  in  the  best  fit  of  the  curves  at  high  and  medium  fre¬ 
quencies.  The  mean  wind  speed  and  wind  direction  measurements  were  performed 
about  fifty  feet  from  the  receiver  building,  hence  the  wake  effects  of  the 
building  could  be  a  contributor  to  the  factor  of  two  discrepancy  in  the 
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wind  speed  producing  the  best  fit  to  the  experimental  spectrum. 


The  two  spectra  fit  closely  in  the  high  and  medium  frequency 
regions.  At  the  low-frequency  outer-scale  breakpoint,  defined  as 
f  ■  AV|_/2irLe,  the  two  curves  begin  to  separate.  In  the  low-frequency 
region,  both  curves  show  a  flattening  trend,  however,  the  experimental 
spectrum  settles  to  a  lower  value  than  the  theoretical  spectrum.  In 
general,  however,  for  this  day  of  moderate  wind  speed  and  small  wind 
direction  variance  the  experimental  and  theoretical  spectra  showed 
good  agreement. 

Figures  5.6a  and  5.6b  are  the  microtemperature  spectra  for  the 
May  25,  1972  mission.  Figures  5.6a  and  b  correspond  to  5.3a  and  b  in 
that  we  have  the  individual  and  averaged  spectra  as  well  as  the  theory. 
Again,  excess  low  frequency  power  is  also  present  in  the  microtempera¬ 
ture  fluctuations  and  this  supports  the  phase-difference  results. 
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Chapter  6 
CONCLUSION 


A  series  of  experiments  to  measure  the  temporal  power  spectral 
density  of  the  optical  phase  coherence  perturbations  across  a  laser 
wavefront  due  to  atmospheric  effects  was  conducted.  The  basic  prin¬ 
ciple,  and  equations  of  the  underlying  theory  were  presented  and  a 
theoretical  phase  difference  power  spectrum  possessing  a  flat  1„„ 
frequency  asymptote  has  been  introduced.  The  measurement  techniques 
the  data  processing  and  analysis  procedures  were  described.  At 
the  present  time  data  from  the  first  series  of  experimental  measure¬ 
ments  has  been  analyzed  and  data  representing  two  typical  conditions 
was  presented. 

Major  emphasis  was  placed  on  analysis  of  averaged  spectra  to 
obtain  adequate  agreement  with  theoretical  predictions  of  the  spectral 
characteristics.  Since  the  confidence  limits  of  the  resultant  averaged 
spectra  were  dependent  upon  the  number  of  spectral  values  averaged 
together.  It  was  necessary  to  uae  lengthy  data  recording  Intervals 
to  produce  accurate  spectral  curves  at  the  low  frequencies  of  Interest. 
For  example,  to  produce  spectral  values  at  2.44  ,  10“3  hertz  with 
ten  degrees  of  freedom  It  was  necessary  to  record  ten,  6.8  „l„„te 
data  record,  or  a  total  of  68  minutes  of  data.  To  maintain  the 
statlonarity  of  the  characteristics  of  this  atmospheric  turbulence 
generated  phenomenon  It  was  therefore  necessary  to  avoid  gross  changes 
in  meteorological  conditions  over  time  periods  of  several  hours.  This 
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restricted  the  measurements  to  days  of  fairl 
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In  the  low-frequency  region  below  this  frequency  the  power  spectral 

curves  were  in  error  by  an  order  of  magnitude,  the  theoretical 

values  being  less  than  the  experimental  ones. 

The  excess  power  in  the  low  frequencies  of  the  phase- 

difference  spectra  as  compared  with  theory  (Figures  5.3  and  5.6) 

was  traceable  to  excess  power  in  those  same  frequencies  in  the 

microthermal  spectra  (Figures  5.3b  and  5.6b).  Since  temperature 

variations  almost  entirely  govern  refractive  index  fluctuations  for 

optical  and  infrared  wavelengths,  we  determined  that  the  theoretical 

basis  for  the  phase-difference  spectrum  was  not  adequate.  That 

basis  is  the  commonly-used  von  Karman  spectrum  described  by 

(3) 

Reinhardt  and  Collins  .  A  new  model  which  has  more  low  frequency 
power  and  which  better  matches  the  microthermal  spectra  is  the  sub¬ 
ject  of  a  separate  report  by  Greenwood  and  Tara-^ano^10^ . 
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APPENDIX  A 


Calculation  of  the  Experimental  Interference  Pattern 

The  section  of  the  Perkin  Elmer  analyzer  which  forms  the  phase 
interference  pattern  is  shown  schematically  in  Fig.  A.l.  In  this 
diagram  the  secondary  Dall-Kirkham,  a  reflective  imager  or  condenser, 
is  replaced  by  its  optical  equivalent,  a  converging  lens  of  identical 
focal  length  and  relative  aperture  number.  The  energy  illuminating 
the  sampling  apertures  is  considered  to  be  phase  coherent  across 
each  individual  aperture  with  a  time-varying  phase  difference,  <J>(t), 
existing  between  the  illumination  at  the  two  apertures. 

To  calculate  the  image  plane  intensity  distribution  it  is 
necessary  to  calculate  separately  the  image  plane  amplitude  distri¬ 
bution  due  to  the  field  illuminating  each  aperture,  add  the  two 
resulting  complex  image  plane  amplitude  distributions,  and  then  com¬ 
pute  the  resultant  image  plane  intensity  distribution  hy  finding  the 
modulus  squared  of  the  sum  of  the  complex  amplitude  distributions. 
These  calculations  will  now  be  performed  to  yield  the  image  plane 
intensity  distribution  I(x,y). 

The  image  plane  amplitude  distribution  due  to  the  illumination 
of  aperture  one  by  the  field  Dje®1  is 

Ej(x,y)  ■  //Djel8le  lkf  ^^^dCdn  (A-i; 

where  x  and  y  are  the  image  plane  rectangular  coordinates 
k  -  2tt/A 

X  ■  optical  wavelength 
f  ■  focal  length  of  the  convergent  lens. 
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The  limits  of  intergration  extend  over  aperture  one  for  (A-l). 

Changing  to  the  primed  rectangular  coordinate  system  (£'  and  n'), 
whose  origin  lies  at  the  center  of  aperture  one  and  integrating  over 
aperture  one  transforms  (A-l)  into  the  form 

Ei  (x*y)  -  e‘lkxa/2f  //Di  ei8i  e‘ikr ‘  >d£'  dn’  (A-2) 

where  a  *  the  separation  of  the  2  apertures 

V  -  5*  a/2 

n'  ■  n. 

Converting  to  cylind“ical  coordinate  systems  in  both  the  aperture 
and  image  planes  changes  (A-2)  to 

ikracos(0)  d/2  r.  is^-ikf- ‘rptcosWcos^)  + 

Ej  (r  ,G)*e - 2f  1  {  Die  6 

p-o  <P*o  (A-3) 

sin(0)  sin ($)  ]  pdpd<^ 

where  r  and  0  are  the  image  plane  cylindrical  coordinates  whose  origin 
is  on  the  optical  axis 

p  and  4*  are  the  aperture  plane  cylindrical  coordinates  whose 
origin  is  at  the  center  of  aperture  1 
■  p  cos^ 
n'  -  p  sin<t> 
x  *  r  cos0 
y  ■  r  ein6 

d  -  diameter  of  apertures  one  and  two. 

By  assuming  that  D!  and  Si,  the  amplitude  and  phase  of  the  incident 
wave,  are  independent  of  p  and  (A-3)  can  be  simplified  to 

E  fr  0)  -  D,el8ie-lkraC08-^  ^  f  e-ikrPf"lc08(Y)PdPdY  <*“*) 
i'  *  '  *  2f  o»o  Y-o 


74 


for  y  «  6  -  <f>. 

The  lnte8ratIo„  th*  of  Y  -ay  be  performed  by  use  of  the 

identity 


27Tjn(z)  -  i"n  }"  eiz  C0S(Y).iny 


dy 


(A-5 ) 


to  yield 


F.,(r,e)  -  2™1.1<S.  -  ^°sW-)f2 
The  identity  /pVj0  (a)da  -  bd,(b)  p.r.itTevOu.tlon  of  the  Integtal 


p.0  pJ»(krpf  1 )dp  (A-6) 


in  (A-5)  as 


Ei(r,0)  -  2-rrD  1e±  (s  1 
which  can  be  written  as 


i(s  _  !«SL£! Si ?j£l  d 

1  2f 


(?) 


Pi 


(krd/2f ) 


(krd/2f ) 


Mx»y)  -7rl^ei(si  '  (d)2 

where  r  -  (x2+y2)  ^2. 


2J , (krd/2f ) 


L  (krd/2f ) 


)_ 


(A-7) 


(A— 8) 


For  K,  =  kd/2f  and  K2  =  ka/f  the  image  plane  amplitude  distribution 
due  to  the  illumination  only  at  aperture  one  is 

M*.y)  «"D,et(si  ~  (i)2  r2Ji(K.lx2+val^lTI 

2  I  (K7fx:+y2]^)  T  (A-9) 

F°r  Ez(X’y)’  the  ima*e  Plane  amplitude  distribution  due  to  the  illumi¬ 
nation  at  only  aperture  two,  a  similar  derivation  is  possible,  the 
only  exception  being  that  the  quantity  £’  is  defined  as  £•  -  £  ~  o/2 

The  result  is  a  sign  change  in  the  argument  of  the  complex  exponential 
phase  term  yielding 


E2(x.y)  “TTl^ei^S2  +  ^2^'  (dj2 


2Ji  (K,  [x2+y 


2+v2-|l/2 


(K,[x2+y 


it *1 

w 


(A-10) 
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The  Image  plana  Intensity  distribution  which  is  defined  as  the  modulus 
squared  of  the  image  plane  amplitude  distribution,  (Ej  +  E2),  is 

Ux.y)  -  [Ej (x,y)  +  E2(x,y)].  [E,(x,y)  +  E2(x,y)]*  (A-li 

or 


I(x,y)  -  EjEj *  +  E2E2*  +  E,E2*  +  £,£,* 
where  *  implies  complex  conjugate. 

Separately  evaluating  the  above  terms  and  for  simplicity  assuming 


that  Dj  -  D2 
I(x,y)  -  M 


we  derive  the  final 


expression 

[1  +  cos  {K2x  +  (Sj 


S2)}] 


(A-13) 


where  M  ■  4irlDl  2 

For  the  arrengement  used  for  this  experiment  the  pertinent  parameters 
were  a  .  0.0381  meters 


d  ■  0.00381  meters 
A  ■  1.06  x  10  5  meters 
f  ■  0.903  meters. 

Based  on  these  parameters  the  half-width  of  the  Airy  function, 

[  ( 2 J j (Kj [x  +y  ]  )}  /  { K j [x2+y 2 ] ^ } ] 2 ,  to  its  first  zero  is  3.12mm. 
The  length  of  one  complete  period  for  the  cosinusoidal  interference 
fringes  is  0.254«.  Therefore,  there  are  twelve  complete  cycles  of 
the  interference  fringes  within  each  half-width  of  the  Airy  function. 
Equation  (A-13)  also  illustrates  the  important  fact  that  (he  phase 
difference  between  the  apertures,  +(t)  -  S,(t)  -  s2(t),  phase  modu¬ 
lates  the  spatial  position  of  the  cosinusoidal  interference  fringes 
in  the  x-direction  of  the  image  plane.  An  optical  phase  difference 
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of  2tt  radians  corresponds  to  a  .254mm  or  one  fringe  period  displace¬ 
ment  of  the  interference  fringes  in  the  x-direction. 


APPENDIX  B 


Reticle  Error  Analysis 

This  analysis  concerns  the  error  in  the  measurement  of  the 
position  of  the  cosinusoidal  interference  fringes  within  the  central 
lobe  of  the  aperture  diffraction  pattern.  This  error  arises  from 
the  use  of  a  circular  chopping  reticle  such  as  described  in  chapter  3 
instead  of  an  infinite  series  of  parallel  square  waves  of  spatial 
frequency  K2  matching  the  spatial  frequency  of  the  interference 
fringes . 

For  this  analysis  the  assumption  is  made  that  the  center  of 
the  interference  pattern  is  located  at  the  proper  radial  distance 
from  the  center  of  the  reticle  such  that  the  arc  length  of  the  first 
radial  wodge  measured  at  that  radius  is  equal  to  one  period  of  spatial 
frequency  K2  *  ka/f 
where  k  =  2rr/A 

A  =  optical  wavelength  of  the  illumination 

a  *  separation  of  interferometer  apertures 

f  ■  focal  length  of  a  converging  lens  used  to  form  the  inter¬ 
ference  pattern. 

For  the  i-th  wedge  or  cycle  there  are  two  errors;  an  error  in 
the  determination  of  the  absolute  location  or  position  of  that  cycle 
denoted  as  Elfi,  and  an  error  in  the  measurement  of  the  length  of 
the  i-th  interval  denoted  as  E2 


The  positional  error  for  each  cycle  of  the  reticle  is  greatest 


at  that  end  of  the  cycle  farthest  from  the  center  of  the  interference 
pattern.  The  calculation  of  the  positional  error  will  be  performed 
for  this  point  to  yield  the  worst  case  positional  error  per  cycle. 

The  positional  error  will  be  expressed  in  percent  of  the  ideal  length 
of  one  cycle  according  to  the  formula 


E1  i  "  {tSi  ‘  V/S!}  x  100  (B-D 

where  Sj  •  arc  length  to  farthest  edge  r t  i-th  cycle  from  inter¬ 
ference  pattern  center  which  by  definition  is  the  ideal 
location  of  the  i-th  cycle 

■  height  at  which  the  radial  line  along  farthest  edge  of 
i-th  cycle  intersects  tne  vertical  axis 

■  ideal  length  of  one  cycle  for  spatial  frequency  K2  and 
also  the  arc  length  of  the  first  cycle. 

From  the  geometry  of  the  problem  the  quantity  may  be 
expressed  as  «  r  tanffi^  which  by  use  of  the  arc  length  relation¬ 
ship,  »  r0^,  becomes 


H  -  —  tan(0.  )  . 

ei 


(iS-2) 


The  quantity  r  is  the  radial  distance  from  the  center  of  the  reticle 
to  the  center  of  the  interference  pattern.  The  quantity  0^  is  the  plane 
angle  on  the  face  of  the  reticle  and  is  defined  as  the  angle  separating 
the  radial  line  to  the  center  of  the  interference  pattern  and  the 
radial  line  to  the  far  side  of  the  i-th  cycle. 


79 


Substituting  (B-2)  into  (B-l)  gives  the  following  expression 

for  the  positional  error 

E1  i  “  {lS  -lSi/eiltan(01)l/s1}  x  WO* 

Ibe  src  length  of  eeoh  Interval  or  cycle  of  the  reticle  is  eguel, 

therefore, 

St  “  i  sl- 


(B—  3) 


(B-4) 


Equation  (B-3)  can  now  be  simplified  to  the  form 


l»i 


{i[l-t  [tan(i6j)  l/tiSj^lH  *  100 * 


(B— 5  ) 


The  interval  length  error  is  computed  as  a  percentage  of  the 
Ideal  interval  length  according  to  the  formula 

E2,i  “  {lsflVVi11/si1  x  10°* 

Equation  (B-2)  can  be  used  to  reduce  (B-5)  to 

*2 1 .  u-us^vms^n  +iisl_lt“,tei_i>,/ls1ei-i1,>  ’  10°- 

Substitution  of  <B-4)  into  (B-7)  gives  the  equation  for  the  interv.l 
length  error  as 

E  t  -  (i-t[tm(ie1)]/!eln+[lt«.ai-iiei»'leil»  *  10°- 

For1  the  particular  esse  of  this  experiment  the  reflective 
reticle  consist,  of  940  co^lete  cycle,  of  sp.tl.l  frequency  K2. 

E.ch  cycle  coneiet.  of  .  one-h.lf  period  length  reflecting  surf.ee  and 
.  on.-h.lf  period  -..orbing  surf.ee.  The  gusntlty  Bj  is  egusl  to 
2F/940  or  0.00668424  redisns.  Initial  attempts  to  calculate  the  err 
b,  use  of  the  power  series  approximation  for  the  tangent  of  angle  <>1 
showed  the  necessity  to  Include  the  fifth  power  of  the  angle.  Sub¬ 
sequently  a  short  computer  program  was  written  to  calculate  both 
errors  by  use  of  the  computer's  tangent  function.  Table  B.l  contains 
a  tabulation  of  both  errors  for  the  first  twenty-five  reticle  cycles. 


(B-6) 


(B-7) 


(B-8) 
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Table  B.l 

Error  Due  to  Circular  Reticle 


Interference 

larcent 

Percent 

Fringe 

Positional 

Interval 

Number 

Error 

Length  Error 

1 

1. 489x10" * 

1.489x10"* 

2 

1.192xl0"2 

1.043x10" 2 

3 

4.022xl0"2 

2.83x10  "2 

4 

9.53xl0"2 

5.51xl0"2 

5 

1.86x10"* 

9. 09x10" 2 

6 

3.22x10"* 

1.36x10" 2 

7 

5. 11x10“ 1 

1.89x10"* 

8 

7.63x10"* 

2.52x10"* 

9 

1.09 

3.24x10"* 

10 

1.49 

4.05x10"* 

11 

1.99 

4.95x10"* 

12 

2.58 

5.94x10"* 

13 

3.28 

7.02x10"* 

14 

4.10 

8.19x10"* 

15 

5.05 

9.46x10“* 

16 

6.13 

1.08 

17 

7.35 

1.23 

18 

8.74 

1.38 

19 

1. 03x10 1 

1.55 

20 

1. 20x10 1 

1.72 

21 

1.39x10* 

1.90 

22 

1.60x10* 

2.09 

23 

1.83x10* 

2.30 

24 

2.08x10* 

2.51 

25 

2.35x10* 

2.73 
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Table  B.l  shows  that,  except  for  the  first  cycle  where  both  errors  are 
equal  and  are  1.489  x  10  3  percent,  the  positional  error  is  greater 
than  the  interval  length  error.  As  explained  in  the  body  of  the  text, 
the  highest  order  fringe  or  cycle  to  effect  the  experimental  measure¬ 
ment  is  the  seventh  cycle.  For  the  seventh  cycle  the  positional  error 
is  0.5113  percent  and  the  interval  length  error  is  0.1894  percent  to 
give  a  worst  case  error  of  0.7007  percent  of  one  interval. 
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of 

Rome  Air  Development  Center 


RADC  is  the  principal  AFSC  organization  charged  with 
planning  and  executing  the  USAF  exploratory  and  advanced 
development  programs  for  electromagnetic  intelligence 
techniques ,  reliability  and  compatibility  techniques  for 
electronic  systems,  electromagnetic  transmission  and 
reception,  ground  based  surveillance ,  ground 
communications ,  information  displays  and  information 
processing .  This  Center  provides  technical  or 
management  assistance  in  support  of  studies,  analyses, 
development  planning  activities ,  acquisition,  test, 
evaluation,  modification,  and  operation  of  aerospace 
systems  and  related  equipment. 


